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Abstract 

The ground state nature of the Falicov-Kimball model with unconstrained hopping of electrons 
is investigated. We solve the eigenvalue problem in a pedagogical manner and give a complete 
account of the ground state energy both as a function of the number of electrons and nuclei 
and as a function of the total number of particles for any value of interaction U £ M. We 
also study the energy gap and show the existence of a phase transition characterized by the 
absence of gap at the half-filled band for U < 0. The model in consideration was proposed and 
solved by Farkasovsky HJ for finite lattices and repulsive on-site interaction U > 0. Contrary 
to his proposal we conveniently scale the hopping matrix to guarantee the existence of the 
thermodynamic limit. We also solve this model with bipartite unconstrained hopping matrices 



in order to compare with the Kennedy-Lieb variational analysis | KL 
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1 Introduction 

The Falicov-Kimball model [FK| of correlated (spinless) electrons on a lattice A C Z d is governed 



by the second-quantized Hamiltonian 



n = - tx,v c l°y + 2 u J2 w * 4 c * (i-i) 

x,y£A xGA 
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where c' x and c x are the creation and annihilation operators for the electrons at the site x, and 
t Xt y is the matrix element for hopping between two sites. Each lattice site may be occupied by at 
most one fixed nucleus which interacts with the mobile electrons via a on-site interaction U. The 
interaction is repulsive when U > and attractive if U < 0. The occupation number w x is or 1 
according to whether the site x is occupied or not by a nucleus. 

Although the w variables are referred as nuclei, this description is only one of many different 
model's interpretations. One may think of the Hamiltonian ( |1 . 1| ) as an approximation to the 
Hubbard model in which one of the two types of electrons are kept static. See references ||KL| , [K| 
for a survey on this. 

Compared to the Hubbard model, the Falicov-Kimball model has the advantage to be reducible 



to a single particle Hamiltonian. Since there are no interactions between the electrons, (IT) is the 
second quantized form of the Schroedinger Hamiltonian on (A): 

H = -T + 2UW (1.2) 

where T is the self-adjoint operator with matrix elements t xy = t yjX and W is the multiplication 
operator by w x , i.e., the diagonal matrix W = diag ({w x } xe \). Denoting L the number of sites in 
A (L = |A|), T and W are L x L matrices. 

To discuss some of the known results on the ground state of (171), let A be the union of two 
disjoint sub-lattices A and B and let us assume that t x>y vanishes when x and y belongs to the 
same sub-lattice. It is also required that A is T-connected in the sense that T is an irreducible 
matrix (it cannot be written as a direct sum). For example, these conditions are met if T is the 
usual nearest neighbor matrix 

t = S 1 if I*- 3/1 = 1 (L3) 
x ' y 1 otherwise 

t > 0, since any AcZ d connected by their links is T-connected. Note that T defines a bipartition 
of A: x = (xi, ■ ■ ■ , Xd) belongs to A or B according to whether x\ + • • • + Xd is an even or odd 
number. 

Now we introduce some notations e definitions. Let 

N e = E 4 C * and N n =Y,W x (T4) 

xeA xeA 

denote the total number of electrons and the total number of nuclei, respectively, and set M = 
N e + N n . For each fixed configuration of nuclei w = {w x } x£ a and N e , 7i has a ground state energy 
E(w, N e ). From this we define two other kinds of ground state energies: 

E^(N n ,N e ) = min E(w,N e ) (1.5) 

w.N n fixed 

and 

EP\M) = min E^(N n , N e ) . (1.6) 

N„,N e :J\[ fixed 

We shall also denote by Wq, Q C A the configuration of nuclei with w x = 1 if x G Q and otherwise. 
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Kennedy and Lieb [|KL|| have proven the following inequality 

E {1) > ~ Tr (T 2 + U 2 ) 1/2 + U M - X - U | A| (1.7) 
and concluded, from this, the following statements about the ground state of (|1 . 1|) : 

A. Let (U,Af) be such that U > and Af > |A|, with |A| the number of sites in A. Then the 

ground state energy has a minimum value which saturates the inequality ( |1.7| ) at exactly 
the boundary Af = |A|. The ground state is doubly degenerated with the minimum of E 
attained at either (w,N e ) = (uu, \B\) or (w,N e ) = (wg, |.4.|). 

B. Let (U,Af) be such that U < and Af < 2\A\. Then the ground state energy E^ has a 

minimum value given ( |1.7| ) at exactly the boundary Af = 2\A\. The minimum of is 
attained at N n = N e = \A\. The ground state is unique unless |^4| = \B\ = |A|/2. 



Remark 1.1 Notice that, if \A\ ^ \B\, the ground state energy has an asymmetric behavior when 
it passes from the A to the B condition These two ground state energies coincide with at 
J\f = |A| only when \A\ = \B\ = |A|/2. In this situation the minimum of E^> is attained at 
N n = N e = |A|/2. The ground state is doubly degenerated at the nuclei configuration and w®. 
The point N e = |A|/2 is called half-filled band because half of the available eigenstates of H are 
filled with electrons. Due to the fact that nucleus and electron have opposite charges, N n = N e is 
often said to be the neutral point of E^ . It is also worth of attention that, when T is given by 
the nuclei occupy alternate sites on the lattice forming a checker-board. We stress that A and 
B results are independent of the hopping matrix entries (irrespective to their signs and strengths) 
provided t x ^ y vanishes if x and y belongs to the same sub-lattice. 



Kennedy and Lieb |KL] have also proven existence of energy gap at half-filled band. They 



propose two gap definitions according to whether (|1.6|) or (|1.5|) is chosen as the ground state 
energy. Let fi(M) := E^(J\f + 1) — E^(J\f) be the chemical potential and let 

d«(A0 := /i(A0 - /i(J\f - 1) = E^\M + 1) + E {1 \AT - 1) - 2£«(A0 (1.8) 

be the chemical potential discontinuity at N '. The system is said to have a gap of the first kind 
at J\f if there exist a constant E\ > 0, uniformly in A, such that > e\. Analogously, the gap is 
said to be of second kind at (N n , N e ) if 

S 2 \N n ,N e ) := E^ 2 \N n} N e + 1) + E®(N n ,N e - 1) - 2E^(N n ,N e ) > e 2 (1.9) 

for some constant e 2 > 0, uniformly in A. 



C. Let T be an irreducible matrix such that ||T|| < r and t X:V > S for all non-vanishing entries, 
uniformly in A, for some constants r < oo and 5 > 0. Under the conditions of item A, 
there is a first kind gap at Af = |A| and second kind gaps at (N n ,N e ) = (\A\,\B\) and at 
(N n , N e ) = (\B\,\A\) with e 2 > S\ > depending only on 5, r, U. 
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D. Let T be as before. Under the conditions of item B, there are first kind gaps at Af = 2 \A\ and 
at Af = 2\B\ and a second kind gaps at (N n ,N e ) = (\A\, \A\) and at (N n ,N e ) = (\B\, \B\) 
with e 2 > E\ > depending only on 5, r, U. 

The Falicov-Kimball model remains object of intense investigation. For a comprehensive survey 



of this model we quote a recent review by Gruber and Macris [|GM| and references therein. Up to 
the present time, very few results about the ground state do not require half-filled band or neutral 
point condition. Other model restrictions which are frequently assumed concern with the spatial 
dimension d and the strength of interaction U. The model which has been investigated in more 
detail has T given by (|1.3| ) with A C Z (d = 1) (see ||GM|| for many intriguing questions on this). 



Also sufficiently large \U\ are assumed very often. There are not many results when T does not 
vanish at same sub-lattice. 

In view of this scenario it seems worth to investigate one particular model where the eigenvalues 
and eigenvectors of H can be computed explicitly for all A. This program can be accomplished 
if the hopping matrix T has all non diagonal elements equal to a constant. We call this hopping 
matrix mean field matrix in analogy to what is usually called mean field in classical spin systems. 
We warn the reader that the Falicov-Kimball model with mean field hopping matrix cannot be 
confused with the mean field approximation of the Hubbard model. 

The model Hamiltonian ( |1 . 1[ ) with mean field hopping was introduced and solved by Farkasovsky 
JF|, who called it simplified Hubbard model with unconstrained hopping of electrons. Its Hamilto- 
nian is given by (|1.1|) or equivalently by ( |1.2|) with the hopping matrix given by (|A| = L) 

2t 

t*,y = T (1.10) 

for all x, y G A with x ^ y and otherwise. 

Farkasovisky's solution of this model is obtained by a canonical transformation U that diago- 
nalizes H . Unfortunately, there is no mention in reference [FJ] on how could U be computed. So, 
one of the purposes of this paper is to present a simple derivation of the eigenvalue problem of H . 

What makes Q1.10 ) solvable is the fact that the energy does not depend on the nuclei configu- 



ration w = {w x } x£ \: 

E{w,N e ) = EW(N n ,N e ) (1.11) 

since the probability of hopping from any site x to any other site y is evenly distributed. In view 
of this fact the model is not suitable to describe the positional order of nuclei induced by electrons 
(cristalization). However many questions about the ground state energy 

£(1) 

function of Af, 

can be formulated and answered. We believe that these answers expresses general features of the 
Falicov-Kimball model. We also believe that the solution of the eigenvalue problem can be useful 
in studying disordered systems. 

Contrary to reference [0 we have made the hopping matrix size dependent (t X}y — > as L —>■ oo) 
in order to guarantee the existence of thermodynamic limit. This simple fact allow us a complete 
description of the ground state for any value of the model parameters U, L and Af (Farkasovsky 
discuss only the repulsive U > case for finite L). 

The following summarizes our analysis. 
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Theorem 1.2 Let 7i be given by 1\ ) with the hopping matrix T given by Then, for 

L > L Q with L sufficiently large depending on U , we have: 

I. IfU> 0, the ground state energy is a monotone increasing function of Af for all < TV < 

2L. The minimum value of E^ is attained at the following points: 

a. {N n , N e ) = (Af - 1, 1) if < Af < p*L with p* = min (1/2 + l/(4u) , 1); 

b. (N n ,N e ) = ([p*L),Af—[p*L]) if[p*L] < Af < L ([x] means the integer part of the number 
x e R); 

c. (N n , N e ) = (L, TV - L) if L < TV < 2L; 

Moreover, E 1 is constant in the domain [p*L] < TV < L (case b). 

II. IfU < 0, the ground state energy is a monotone decreasing function of N for all < TV < 

2L. Except by 1/L corrections, the minimum value of E^ is attained at the neutral point 
(N n ,N e ) = (TV/2, TV/2) 



Remark 1.3 Monotone behavior of the ground state energy 

£(1) 

seems, from our analysis, typical 

of the Falicov-Kimball model. This property would explain the minima tendency to be attained at 
the domain boundary in A and B results. The derivation of this property for nearest neighbor 
hopping matrix in d = 1 dimension is subject of a forthcoming paper \ML\j . 



Theorem 1.4 Under the condition of Theorem EI, the mean-field Falicov-Kimball model has a 



first kind gap S 1 ' = 2U at the half-filled band Af = L, if U > but no gap of first kind occurs if 
U < 0. 

In addition, there are three gaps of second kind at (N n , 1), (N n , L — N n ) and (N n , L — N n + 1) 
for all nuclei number 0<N n <LifU >0 and at (N n , 1), (N n , N n ) and (N n , N n + 1), < N n < L, 
if U < 0. The gap values are listed in and (fOJ). 



Remark 1.5 Absence of energy gap when U < does not contradict Kennedy-Lieb's result. We 
note that the elements of the mean-field hopping matrix are not uniform in A. We also observe 
that this hopping matrix violates the bipartition condition required in item C. As we shall see, this 
condition does not affect our results on the gap which clearly indicate the existence of a phase 
transition. 



Finally, we illustrate the Kennedy-Lieb results state above with an exactly solvable model. Our 
lattice A may be arbitrarily decomposed into two sub-lattices A and B. The model, called by us 
bipartite mean field, has non-vanishing hopping matrix elements t XjV equal to a constant if x, y 
connects the two sub-lattices. The Hamiltonian is thus given by ( |1 . 1| ) with 

f ^ if x, y is not in the same sub — lattice , . 

x ' v I otherwise 
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The results of itens A D may be appreciated in light of the eigenvalue problem for the present 
model. In particular, one sees that the equation (jl.71 ) becomes equality if one of the the sub-lattices, 
let say A, is constrained to be filled with nuclei (|*4| = N n ) for < TV < 2L when U < and for 
\M/L-l\ > 1/2 when U > 0. 

This paper is organized as follows. In Section |2| the eigenvalue problem of H with mean-field 
hopping matrix ( |1.10 ) is solved. Sections |] and [| give a complete description of the ground state 
energy and energy gaps for this particular model. This analysis is repeated for the bipartite mean- 
field model in Section [|. We finally draw some conclusions in the Section ||. Appendix [A| reviews 
some basic features about circulant matrices used in this text. 



2 The Eigenvalue Problem 

This section is devoted to the eigenvalue problem of the Schroedinger operator ( |1.2j ) with mean 
field hopping matrix ( |1.10|) . 

For convenience, we write the interaction constant U = t u and factorize 2 1 from the Hamilto- 
nian matrix H. This factor is just a scale to the energy and may be fixed equal to 1 without any 
loss of generality. Equation (|1.2|) reads 

H = -T + uW . 

H = H(w) is a L x L matrix with matrix elements h X}V = — 1/L if x ^ y and h X)X = uw x which 
depends on the nuclei configuration w. 

Now we claim that the spectrum of H given by the zeros of the characteristic polynomial 

P(A) = det (H - X I) (2.1) 

is independent of the nuclei configuration. 

To see this, let w and w' be two nucleus configurations differing by a permutation ir, i.e., 
w' x = w n ( x y If II denotes the permutation matrix determined by n, we have II -1 Til = T and 

H(w') = W 1 H(w)U. (2.2) 

Note that the nearest neighbor matrix T given by ( |1.3|) is not invariant by the similarity transfor- 
mation (|2.2|) . 

As a consequence of this fact det (H(w') — A I) = det {H(w) — A I) and the eigenvalues of H 
depend only on the total number of nuclei N n . We shall drop the subscript n of N n in the sequel 
of this section. 

To compute Q2.1D let us pick the simplest matrix H. For any positive integers L and N, L > N, 
let J i , J2 and K be given by the following matrices 

J\ = (u + l/L) I N -{1/L)1 N 

J 2 = {1/L)Il- N -{1/L)1l- N (2.3) 
K = -(1/L)1 N , L _ N 



6 



where Ir is the identity matrix of size R, 1 r,s is the R x S matrix with all elements equal to 1 
and 1 R = M RtR . We set 

* = < 2 - 4 > 

Note that H = H(wq) with the nuclei configuration u> = (1, 1, . . . , 1, 0, . . . , 0). 

In view of the above, the spectrum of H(w) can be read from the spectrum of H given as 
follows: 

Proposition 2.1 The spectrum of H consists of the set {Ai, A2, A + , A_} where 

Ax = l/L 

A 2 = l/L + u (2.5) 
A± = l/L + [u-l±A(N/L)]/2 

with 



A(p) :=^(u + iy-4up } (2.6) 

have multiplicities L — N — 1, N — 1 and 1, respectively. 

The corresponding eigenvectors x = (u, v) £ C w x C L ~ N are Xij = (ujv^ , 0) , j = 1, . . . , iV — 1, 
x 2,fc = (0 , u L _ N ^ k ) , k — 1, • • • , L—N — l, and x± = (1 , u ± ) where = (0, . . . , 0) is the zero-vector, 
1 = (1, . . . , 1) is the one-vector, 

u + l- 2N/L ± A(N/L) 
U± " 2(1 - N/L) 1 (2 ' 7) 

and Urj = (1, u R , io R \ . . . , u R !:j ),j = 1, ... R — 1, with ur = e 2m ^ R . 



The proof of Proposition |2.1| will follow from two lemmas. The first reduces the problem to 



matrices of the form a I + (3 1 . The second describes the properties of such matrices. 
Lemma 2.2 

det (H - XI) = det (J x - XI) ■ det [( J 2 - XI) - K T (J t - XI)' 1 k] . (2.8) 

Proof. Let 



M 



A K 
K T B 



be such that A, B and K are N xN, (L — N) x (L — N) and x (L — N) matrices, respectively, with 
A symmetric and invertible. M can be brought to a block diagonal form by using a non-unitary 
transformation 

( B - K T A~ X K ) = ( -K T A- 1 1 ) ( K T B ) ( V)' (2 " 9) 
Taking the determinant of both sides of equation (|2.9|) with A = J\ — XI and B = J 2 — XI gives 



2.2 . □ 
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Lemma 2.3 1. The linear space C of all R x R matrices of the form a I + (3 1 is closed under 
sum, product and inverse operations. 



2. The inverse of a matrix C = a I + j3 1 is given by 



c- 1 







1 



a a(a + j3R) 
provided a/0 and a ^ —(3R. 
3. Any matrix C = a I + (3 1 can be reduced to a diagonal form 

F' 1 CF = diag {a + /3R, a,..., a} 
where F is the Fourier matrix 



(2.10) 



Fk,i '■- 



R 



(fc-i)(i-i) 

R 



k,l = l,...R 



with uj r = e 27Tt/R . 

In other words, C has a simple eigenvalue Ai = a + (3N corresponding to the eigenvector 
1 = (1, 1, . . . , 1) and a (R — \)-fold eigenvalue A2 = a associated to eigenvectors urj = 



4. detC = {a + f3R)a 



R-l 



The proof of Lemma 2]B is in the Appendix 

Proof of Proposition \2.1\ We use Lemma [2.2| to reduce the calculation of the characteristic function 
-Po(A) = det(H — XI) to algebraic calculations of matrices of C. 

we have 



Using item 1. and 2. of Lemma |2.3| and definitions 

1 



(Ji - A/)" 1 



u+l/L-X N ' (u + l/L- X)(u- (N-l)/L- A) 
and, in view of 1 RtS 1 s ,r = S 1 R , 

N/L 2 



1/L 



1 



N 



RT^-X^K 



u-(N -1)/L- X 



■ L-N 



which gives 



( J 2 - XI) - K T (J, - Xiy 1 K = (1/L - A) I L ^ N - (M t.y L lW A r )/jL > *l-n 

u — (J\ — 1) / L — X 



:2.11) 



We then use item 4- of Lemma EO and equations (|2 .3|) and (|2.11|) to write 



det(Ji - XI) = (u-{N- l)j L - A) (u+ 1/L - X) 



N-l 



[2.12) 
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and, with a little of algebra, 



det ( (J 2 - XI) - K 1 ( Ji - XI)' 1 K 



l/L-A) 



L-N-l 



u- (N -1)/L- A 



(2.13) 



where 



i?(A) = A 2 -(M-l-2/L)A-n(L-A^-l)/L-l/L+l/L 2 
= (A-A+)(A-A_). 



The eigenvalues of H then follows from (|2.12|) , (|2.13| ) and Lemma ^L2 



The corresponding eigenvectors of Hq can be obtained from item 3. of Lemma |2.3| . A vector 
x in C L is written as (u, v) where u e and v G C L ~ N . Using Lemma 22 and the notation 
of Proposition |2.1| , the eigenvectors corresponding to the Ai and A 2 are, by inspection, given by 
(u/vj, 0), j = 1, . . . , N — 1 and (0, u L _ N ^ k ), k = 1, . . . , L — N — 1, respectively Notice K u N j = 
for all j = 1, . . . , N — 1 in view of the property ( |A.7|) in the Appendix [A|. By inspection also, one 
can verify that the eigenvectors corresponding to A± are (l,u±) with u± given by ( [2.7| ). 

This concludes the proof of Proposition |2.1| . □ 



3 Ground States 

For fixed nuclei configuration w, let fix < ^ < • • • < A*i be the sequence of ordered eigenvalues of 



1.2), counting multiplicities, and let N e be the total number of electrons in the system. 

Due to the exclusion principle, each energy level (eigenvalues fij) can be occupied by at most one 
electron. In the ground state the first N e lowest energy levels of H are occupied by electrons and 
the L — N e remaining levels are left empty. The ground state energy E of the model Hamiltonian 
( |1.1| ) at (w, N e ) is thus determined by the sum of the lowest eigenvalues: 

N e 

E(w,N e ) . (3.1) 



We shall here prove Theorem [L2| on the ground state energy of the Hamiltonian ( |1 . 1|) with 



mean field hopping matrix ( |1.10| ) whose eigenvalues have been studied in the previous section. 

Since we have proven the spectrum of H, a{H) = {Ai, A 2 , A + , A_}, depends only on the total 
number of nuclei, we have 

EM(N n ,N e ) =E(w,N e ) 

for all w. Consequently, we shall only be concerned with the minimum value of E^ for fixed 
number of particles Af = N n + N e . 

Our discussion about ground state energy begins by the eigenvalues ordering. We observe that 
the spectrum a(H) depends on the nuclei density p n = N n /L through the function 

A 2 (p) := (u+ l) 2 -Aup (3.2) 
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which, in view of < p < 1, satisfies 

(M-i) 2 <A 2 <(M + i) 2 . 

Since A± = 1/L+ (u — 1 ± A)/2, this implies 

1/L + (w - 1 + | |u| - 1 1) /2 < A+ < 1/L + (w + |u|) /2 

and 

A_ < 1/L + (u- 1 - | \u\ - 1 1) /2 
from which the eigenvalues ordering follows. 
We distinguish two cases: 

I. if u > we have 

A- < Ai < A + < A 2 ; 

II. if u < we have 

A- < A 2 < A+ < Ai . 



Proof of part I of Theorem [7j]. Let u > and M < L. In this case we have 

Ml = A_ , /i2 = • ■ • = fJ'L—Nn = Ai , flL-N n +l = A + , fiL-N„+2 = ' ' ' = = A 2 (3.3) 

and the ground state energy is thus given by 

£( 2 )(iV n ,iV e ) = A_ + (JV e -l)Ai 

= A_ + (AT-iV n -l)Ai (3.4) 
= M/L + g( Pn ) 

where 

:= i (« - 1 - A(p)) - p (3.5) 

and A as in ( |3.2|) . We remind that, under the condition Af < L, the multiplicity L — N n — 1 of Ai 
is greater than or equal to the number of electrons N e — 1. This explain why we have taken only 
two eigenvalues in ( p.4|) . 

Equation ( ^4|) and (|3~5|) reduces the problem of minimizing E^(N n ,N e ) = E^ 2) (N ni M - N n ) 
with respect to N n to the problem of minimizing g as a function of p n . The function g> : [0, 1] — > 1R 
has the following properties: 

1. g(p) < 9(0) = -1 with 0(1) = -1 if u > 1; 

2. attains to a unique minimum value at 



(note that p* = 1 if u < 1/2); 



P =min( 2 + 4^' 1 ' 
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3. the minimum value is given by 



, ,* / u-2 if < it < 1/2 , , 

^ ) = \ -l-l/(4r*) if „>l/2, (3 - 6) 

Items 1, 2 and 3 can be obtained explicitly. Note that the derivative 

M = W) ~l (3.7) 

vanishes only at p* provided u > 1/2. If it < 1/2, <? attains its minimum value at the domain 
boundary p = 1. 

From equation ( |3.4j ) and properties 1,2 and 3 we conclude that the ground state E^(Af) 
has only one eigenvalue A_ occupied by one electron if Af < p* L and there will be more electrons 
contributing to the ground state if Af > p* L. Note that A_ = A_(p„) is negative if p n < p* provided 
L > 2. This follows from the fact that A_(p n ) < A_(p*) = 1/L+ (u - 1 - A(p*))/2 = l/L - 1/2 
since, from ( |3.7|) , A(p*) = u. We shall always assume L sufficiently large. 

In summary, 

F (i)fAA_J M(JV-1)/L) if AT<p*L 

l7VJ ~ \ Af/L + g([p*L]/L) if p*L<Af<L {6 *> 

where [x] means the closest integer of x £ R. We emphasize that E^'(Af) is a monotone increasing 
function of Af, since A_ is a monotone increasing function of the nuclei number, and piecewise 
constant if p* L < Af < L. 

Now, let m > and Af > L. The ground state energy E^(N n , N e ) in this case is given by 

E^(N n} N e ) = \_ + (L-N n -l)X 1 + X + + (N e ~(L-N n + l))X 2 

= A_ + (L-iV n -l)A 1 + A+ + (A/'-L-l)A 2 (3.9) 
= (Af-N n )/L+(Af-L)u-l 

which attains to a minimum value at (N n , N e ) = (L, Af — L) corresponding to the maximal nuclear 
occupation. We thus have 

(Af) = {Af - L) (l/L + u) - 1 (3. 10) 

for L < Af < 2L. The combination of (|3.8f) and ( p,10|) leads to a monotonic increasing ground state 
energy as a function of Af in the whole domain. 

This concludes the analysis of the case I. □ 



Proof of part II of Theorem l.i. Here, we have 



V\ — A- , P2 — ■ ■ ■ — ^N n — A2 , PV„+1 — ^+ ) V>N n +2 ~ " " " — ~ A 2 ■ (3-H) 

We shall divide our analysis in two sub-cases. Let u < and N e < N n . In view of the eigenvalue 
ordering, we have 

E^(N n ,N e ) = A_ + (iV e -l)A 2 

= A_ + (JV- iV n -l)A 2 (3.12) 
= Af/L + u(Af-l) + h(p n ) 
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where 

h(p)=g(p)-uLp. (3.13) 

Once again, the problem of minimizing E^ 2 \N n , N e ) under the condition N e < N n is reduced 
to the problem of minimizing h in the domain 1/L < p < 1. Note that the lower limit 1/L is due 
to the fact the equation ( |3.9| ) is valid only if 1 < N e < p n L. The following description of h gives 
the answer to this problem. 

Proposition 3.1 Given u < 0, there is Lq = Lq(u) such that, for all L > Lq, h : [1/L, 1] — > R 



given by ( \3.13\ ) is a strictly monotone increasing function of p. 



Proof. Proposition |3J] is implied if h'(p) > for all p in [1/L, 1]. From ( [3.13j ) and ( |3.7|) this is 
implied by 

u 

m >uL+l. (3.14) 

To show that condition (|3.14|) holds for any value u G (— oo, 0), we must break this open interval 
into three pieces: (-oo, -1 - 1/L) U (-1 - 1/L, -1 + 1/L) U (-1 + 1/L, -1/L). Note that we have 
to take L large in order get u arbitrarily close to 0. This leads to L > — 1/u. 

Let us assume —1 + 1/ L < u < —1/L. Using the inequality A(p) > u + 1, equation ( E>.14| ) can 
be replaced by a more restrictive condition 

u 

> uL + 1 . 

u + 1 

This is equivalent to Lu 2 + Lu + 1 = L(u + 1 — 1/L)(u + 1/L) < which is always true for 
u G (—1 + 1/L, —1/L). 

Suppose now u < — 1 — 1/L. Then we analogously use A(p) > — (u + 1) to replace ( |3.14| ) by a 
more restrictive condition 

Lu 2 + (L + 2)u + l > 0. (3.15) 

Since the roots u = —1/2 — 1/L ± [1/4 + 1/L 2 ] 1 / 2 of the quadratic equation are close to —1 and 0, 
d3T| ) is always satisfied if u < -1/2 - 1/L - [1/4 + 1/L 2 ] 1/2 <-l-l/L. 

It remains to verify condition ( |3.14| ) for — 1 — 1/L < u < — 1 + 1/L. In this case we use p > 1/L 



to replace (|3.14j ) by 

>uL + l. (3.16) 



A(l/L) ^( M + 1)2 _4uAL 
Since u = 0(1) and (u + 1) = O (1/L), the left hand side of this inequality is O (Vl) and the right 



hand side O (L). Recall u < —1/L and both sides of ( |3.16| ) are negative. Therefore, we can always 
find L sufficiently large so that condition (|3.14j ) is verified concluding the proof of the Proposition. 

□ 



Proposition 3.1 implies that the ground state energy under the condition iV e < N n has a 



minimum value at N n = N e = Af/2 given by 

E^(Af) = Af/L + u(Af - 1) + h(M/(2L)) (3.17) 
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Finally, let us consider u < and N e > N n . The ground state energy in this case is given by 

E&(N n ,N e ) = X^ + (N n -l)X 2 + X + + (N e -N n -l)X 1 

= Af/L - 1 + (Lu - l)N n /L {6 8J 

Since (N n , Af—N n ) is a monotone decreasing function of p n = N n / L, it attains to a minimum 
value at N n = N e — 1 = Af — N n — 1 given by 

E^(Af) = Af/L -l + (Lu-l)(Af- 1)/(2L) . (3.19) 



The ground state for the case II is thus given by 



EW(Af) = min (E^\Af) , E^\Sf) 



u 



Af m (3.20) 



2^-^ + 2L 2 

where m := max (2 - 1/L, 1 + A(jV/(2L))). 

Note that -E^ 1 -* is a monotone decreasing function of Af if L is sufficiently large. We point out 
that the minimum of E^> is attained at approximately neutral (symmetric) point N n + 0(1/ L) = 
N e + 0(1/ L) = Af/2 for all values of u < and < Af < 2L. 



This concludes the proof of Theorem 1.2. □ 



4 Gaps 



Let {ei(p)}i>o be the sequence of ground state energy density functions defined by 

/ ^ E^(pL) 
edp) = z 



(4.1) 



and let e(p) = lim^oo e^(p) be its thermodynamic limit. 

From our previous analysis, the ground state energy E^ can be written in the form (see 
equations Q), (PHQ and (Q)) 

E {1) (pL) = a (u, p, L) + at(u, p) + a 2 (u, p) L (4.2) 

where, for fixed u, ao, a\ and a 2 are bounded functions of p G [0, 2] with ao — > as L — >• oo. 
Hence, given two positive integer numbers L, M, L < M, 



E^(pL) E^(pM) 



M 



< 



a (L) a (M) 



L 



M 



+ ai 



1 1 
L ~~ M 



< C- 



L 



(4.3) 



for some constant C = C(u) < oo, uniformly on L, and this shows that {ei}L>o forms a Cauchy 
sequence which converges to e(p). 

It also follows from O, (|3TT0| ) and (|3T20D that: 
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a. if u > 0, we have 



e(p) = { ° if °" P - 1 (44) 

6[P) \ {p - l)u if 1< p < 2 . l4 ' 4j 



( Note that a 2 = when < p < 1); 
b. if m < 0, we have e(p) = u p/2 for all < p < 2. 

By a similar computation, one can show that the sequence of chemical potentials {pl(p)}l>o 
given by 

MP) := £<»(pL + 1) " ^\ P L) = ^ + V^M (4.5) 

converges to the derivative of the ground state energy density e'(p) for all p G [0, 1) fl (1, 2] provided 
u > 0. As a consequence, there exists a positive gap of first kind at p — 1: 

lim p(p) - lini /i(p) = u . (4.6) 

For u < the sequence of chemical potentials {pl(p)}l>0) converges to e' L (p) = u/2 for all 
p G [0,2] and there is no first kind gap in this regime. 

Now, let us consider the second kind gap. From the definitions ( |1.9| ) and ( |3.1| ) we have 

= Pn e +i - PN e (4.7) 

(For this, recall (|TTlD ). 

For u > 0, it follows from the equation (|3.3|) several cases: 

1. ^(iV n ,l) : ~ U+1 + A ^ 

2. ^(iV n ,iV e ) =0 " if KN e <L-N n 

3. d( 2 )(iV n ,L-iV n ) = M ~ 1 t A(Pn) (4-8) 

4. rf( 2 )(iV n ,L-iV n + l) : 

5. S 2 \N n ,N e ) =0 " if A^ e >L-iV n + l 

We thus find three gaps for u > 0: second kind gaps occur at N e = 1, L — N n and L — N n + 1 
for each < N n < L. 

For u < equation ( |3.11D gives the following: 

1. ^,i) = u -±1±^eA 

2. rf (2) (iV„,iV e ) =0 " if KN e <N n 

3. ^(iV n ,iV n ) = ~ M ~ 1 2 +A(P " ) (4.9) 

4. ^(iV„,iV n + l) = - M + 1 - A( ^ 

5. c/( 2 )(iV n ,iV e ) =0 " if N e >N n + l 





2 







It - 1 


+ A(p n ) 


u+ 1 


2 

- A(Pn) 




2 
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From these we find that, for u < 0, second kind gaps occur at N e = 1, N n and N n + 1 for each 



< iV n < L and this concludes the proof of Theorem |LJ. □ 



5 Bipartite Mean Field Model 

We shall briefly in this section analyze the ground state of the Hamiltonian 7i with bipartite 
mean-field hopping matrix T given by ( |1.12| ). 

We begin by solving the eigenvalue problem. Let us decompose the vector space C L , L = |A|, 
into two subspace and C M corresponding to each sub-lattice A and B where N = \A\ and 
M = \B\. The nuclei occupation number in a sub-lattice X C A will be denoted by nx- We shall 
frequently write n = and m = ng. From these notations we have L = N + M and N n = n + m. 

The Hamiltonian H = H(w), acting on C L , can thus be written as a block matrix 

(51) 

where Jx = u dia.g{w X) x G Af}, X = A,B and = — (1/L) £n,m- 

As before, we can find a permutation matrix II such that ifo — n _1 ifn has the diagonal 



block matrices of the form diag{w, . . . , u, 0, . . . , 0}. Hence, in view of Lemmas 2/2 and 2.3 , the 
characteristic polynomial P(X) of Ho is given by 

P(A) = det (#„ - A/) = det A det C (5.2) 

where A = (J4 — XI) and C = (Jg — A/) — fT T (J4 — A/) - if, can be written as 

C=(££) (5.3) 

with 

E = (u - A) I m - a l m 

F = -XI M - m ~ a 1 M _ m (5.4) 
G = —al m M-m 



and 



1 (u — \) N — un 
a= V \{u-\) ■ (5 ' 5) 



We now repeat the calculation given in the proof of Proposition j2.1| . Once again, by Lemmas 



0]and|2J 

det C = det E det H (5.6) 

where 

H = F-G- l E- l G 

a(u-\) (5.7) 

— —A lM-m ^ R M-m 

u — A — am 
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Putting all pieces together, by using item 4 of Lemma [2.3| we have 



detA = (u - A) n (-X) N ~ n 

detE = (u- X)™- 1 (u- A- am) 



detH = -(-A) 



M-m-i X ( u ~ A ) + a H M -m)-XM] 



(5. 



u — \ — am 

The characteristic polynomial -P(A) = detA det E detH can thus be written as 

P(A) = ( u - A)^- 2 (-\) L - N »- 2 5(A) (5.9) 
(recall L = N + M and iV n = n + m) with 

5(A) = A (w - A) - ( — A — u ) ( —A F — u ) . (5.10) 

We have proven the following 



L L J \ L L 



Proposition 5.1 The spectrum of H given by (5.1) is the set {Ai, A2, A +j+ , A +i _, A_. + , A-^} wh 



wre 



Ai = 0, A2 = u and A±,±, the four solutions of the equation 5(A) = 0, have multiplicities L — N n —2, 
N n — 2 and 1, respectively. 



Remark 5.2 Although the eigenvalues X±,± can be evaluated analytically, their expressions are too 
large to be useful. Numerical computation indicates that for any values of parameters p := N/L, 
Va := n l L and r/g := m/L in the interval (0, 1/2), we have the following ordering 

A- _ < Ai < A+_ < A_ i+ < A 2 < A +>+ if u > , . 

A_ _ < A 2 < A+ _ < A_ i+ < Ai < A +>+ if u < V ■ ) 

The solutions of 5(A) = can be easily found if two roots are equal to and u. This equation 
can be easily solved also when the second term of the right hand side of (|5.10| ) is a square of a 
linear function of A times a constant. We shall consider three choices of parameters having these 
properties: 

Example 5.3 1. If p < 1/2 (N < M), rj A = p = N n /L (n = N) and r] B = (m = 0), the 

characteristic polynomial 

S(\) = A (A - u) (A (A - u) + p(l - p)) 
has the following roots: A +i _ = Ai = 0, A_ i+ = A 2 = u, A_ t _ = A_ and A + + = A + where 



u / /w x 2 



X±= 2 ± \\V (5 - 12) 
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2. If p = 1/2 (N = M = L/2) and r] A = r] B = 77 = N n /(2L) (n = m), we have 

S{\) = A 2 (A - uf - (A/2 - (1/2 - r])u) 2 

which gives 

A+,± = (u+l/2±A+)/2 
A_ i± = (u-l/2±A_)/2 

where 



(5.13) 



A ± (r / ) = yJ(u T 1/2) 2 ±4m^ (5.14) 
5. Ifp = N/L, r] A = p/2 (n = N/2) and rj B = (1 - p)/2 (m = M/2), we have 

S(X) = A 2 (A - uf - p(l - p) (A - m/2) 2 
which has the following roots 

A±,± = (5.15) 

where Q± = Q±(p) is given by 

&l=u 2 + 2p(l - p) ± 2^(1 - p) 2 + w 2 p(l - p) (5.16) 
(Note that here N n = n + m = L/2) 

For comparison with the variational method used by Lieb-Kennedy and with the results stated 
in Section [| we shall compute the ground state energy with these eigenvalues. Our analysis here, 
opposed to the previous sections, will be brief and very selective. 

Let N, M, n and m be as in the Example \5. In this case we have A_ < Ai < A2 < A + if u > 
and A_ < A 2 < Ai < A + if m < 0. 

Proceeding as in the Section [| the ground state energy at (N n , N e ) is given by 
I. For u > 0, 

a. if Af = N n + N e < L, we have 



E^(N n ,N e ) = \_(p) + (N e -l)\ 1 
= A_(p) 



(5.17) 



where p = N/L = N n /L; 
b. if M > L, we have 



E^(N n ,N e ) = \_( p ) + (N e -L + N n )\ 2 
= A_(p) + (Af -L)u 



(5.18) 
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Now, the ground state will be minimized with respect to (N n , N e ) with Af = N n + N e 
fixed. Note that the eigenvalue A_, given by ( |5.12| ), is a monotone decreasing function of p 
in < p < 1/2 and a monotone increasing function inl/2<p<l for all «6l This gives 

c. if Af < L, 

E^{Af) = A_(p*) (5.19) 

where p* = min((jV- 1)/L, 1/2); 

d. if N > L 

E^{Af) = A_05) + (Af - L) u (5.20) 
where p = max (1/2, Af /L — 1). 

The minimum is attained at (N n , N e ) = (Af - 1, 1) if TV < L/2, (N n , N e ) = (L/2, Af - L/2) 
if L/2 < Af < 3L/2 and at (N n , N e ) — (Af — L, L) if Af > 3L/2. As the number of sites at 
the sub-lattice A coincide with the number of nucleus, |^4| = N = N n , the minimum value 
of E^ is attained at the partition of A into two evenly divided sub-lattices, \A\ = \B\, for 
L/2 < Af < 3L/2. We recall that the sub-lattice A is completely occupied by nuclei and the 
sub-lattice B is empty in this example. 

II. For u < 0, we have 



a. if N e < N n , 

E^(N n ,N e ) = A_(p) + (N e -l)u (5.21) 



b. if N e > N n , 



E {2 \N n , N e ) = A_(p) + (N n 1) it (5.22) 
Since the minimum of E® occurs at N n = N e = Af/2, this gives 

E^(Af) = -\ + (Af/(2L))+Af/2u (5.23) 



□ 



Let N, M, n and m be as in the Example [57j.[|. Using the fact that, for u > 0, 

\u- 1/2| < A± < u + 1/2 (5.24) 

and Proposition [□], we have p ± = A_ _, (jl 2 = • - - = Pl-n u -i = 0, Pl-n u = A+ _, Pl-n u +i = A_ +, 
PL-N n +2 = ■ ■■ = Pl-i = u and p L = A +)+ . 

We also have 

\u + 1/2| < A± < 1/2 - u , 



if m < 0, which gives /ii = A_ _, p 2 
A*l-i — and pl = A ++ . 



+2 
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Observe that A +i _ < A_ i+ for any «el and < rj < 1/2 because 

A+ + A_ > {2{u 2 + 1/4) + 2\u 2 - 1/4|} 1/2 > 1 . 

Proceeding as in the Section |], the ground state energy E^ at (N n , N e ) is given by 

I. For u > 0, we have 

a. if Af = N n + N e < L, with N e ^ L, 



E^(N n ,N e ) = A_,_(77) + (N e - l)\t 
= ( W -l/2-A_(7j))/2 



(5.25) 



(recall r? = N n /(2L)); 

b. 

c. if Af = N n + N e > L, 



E^(0, L) = A_,_(0) + A+,_(0) = -\ ; (5.26) 



E^(N n ,N e ) = \_ t _(rj)+\ + ,-(Tj)+\- t+ (ri) + (N e -L + N n -i)u 
= (u-1/2- A+(r)))/2 + (N - L)u 



(5.27) 



We note, from definition Q5.14| ), that A CT is a monotone increasing function of rj if w > and 
cr = + ( or if -u < and a = —). On the other hand, A CT is a monotone decreasing function 
of 77 if u > and a = — (or if w < and a = +). Minimizing i?( 2 ) with respect to (iV n , N e ) 
with Af = N n + N e fixed, gives 

d. 

E^(Af) = {u- 1/2 - A_(0))/2 = -1/2 (5.28) 

if Af < L and 

E^{Af) = (u - 1/2 -A + (l/2))/2+ (AT -L)u 

= -1/2 + (Af -L)u 15 J) 

ii Af > L. 

Note that the minimum is attained at (N n ,N e ) = (0,7V) if TV < L and at (N n ,N e ) = 
(L,Af — L) if Af > L. In other words, the ground state has no nucleus occupying the lattice 
sites of A if Af < L and A is fully occupied by nuclei if Af > L. 

II. For u < 0, we have 
a. if N e < N n , 

E&(N n ,N e ) = \_^ v ) + (N e -l)u 

= -{u + 1/2 + A_{rj))/2 + {AT- N n )u { °- 6U> 
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b. if N P > N, 



n ■ 



E^(N n ,N e ) = X^4 V ) + (N n -2)u + X +t 4 V )+X^ + ( V ) 

= -(u + l/2 + A + (r)))/2 + N n u { °- 6i> 

To minimize in this case one has to deal with two contributions to the ground state energy. 
Note in the equations (|5.30|) and ( 5.31|) that there is one term depending on 77 = N n / (2L) 



and other depending on N n . As N n varies, these two contributions moves the ground state 
energy into opposite directions. Despite this, we claim that, if L is large enough, the second 
contribution always dominates and the following description holds (see Section ^| for similar 
analysis) : 

c. 

£ (1) (A0 =-~-~ max( A + (C), A_(C) ) + (5-32) 

where ( = N / (4L). Note that the minimum is attained at the neutral point N n = N e = 
N/2. In addition, note that A+(C) = A_(C) at C = 1/4 {H = L). Since A CT (C) increases 
or decreases according to whether a is positive or negative, 

E^\U) = -\-\k{Q + ^^u (5.33) 
with o = - if C < 1/4 and a = + if ( > 1/4. 

□ 



Let N, M, n and m be as in the Example \5. 5| .|3|- Note that, for all u £ M different from zero, 
CT (p), given by (|5.16| ), is a monotone function of p, < p < 1/2, which increases if a = + and 
decreases if a — — . This follows from the fact that G± > and 

.-2p 



e±(p) 



l± 20 (2il 



2p) + u' 



where 9 = (p 2 (l - p) 2 + u 2 p{\ - p)) 1/2 . Note, in addition, that 6±(0) = u and e' ± (0) = ±00. 

This behavior of G± allows us to deduce the expressions of the ground state energies E^ and 
E( 2 ) from those already obtained for the Example 0by making the following replacement: 



±1/2 - A± 
±1/2 + A± 



For the item II. c , e.g., 



E^(Af) 



max(9 + (p),e_(p)) 



u 



(5.34) 



(5.35) 



holds for all u < 0, where p = N/ (2L). We recall that the minimum of E^ is attained at N n = N e . 
Since the nuclei number N n = L/2, we conclude that M = N n + N e = L is the only accessible point 
in this example. □ 
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6 Conclusions 



The nature of the ground state of the Falicov-Kimball model with long range hopping matrices 
(mean-field matrices) was investigated. The eigenvalue problem of this model has been entirely 
solved allowing us to obtain the profile of the ground state energy. We have investigated the 
ground state energy 

£(2) 

as a function of the number of nuclei and electrons (N n , N e ), and as 
a function of the total number of particles Af = N n + N e , for all values of the interaction U . This 
analysis revels general features of the Falicov-Kimball model which may be conjecture to be true 
irrespective to whether bipartition and half-filled band hold. The aim of this section is to point 
out some of these features. 

Kennedy-Lieb's results on the Falicov-Kimball model (items A-D in the introduction) require 
the lattice to be bipartite. Let us compute, for comparison, equation ( |1 . 7| ) for the bipartite mean- 
field model. In this case, we have T 2 = (M/L 2 ) In® (N/ 1 L 2 ) 1m which gives (with U = tu and 
2t = 1) 

TW^ + fUe(^, + ^ f ) 



and, by using Lemma 2.3 



(6.1) 



where 

\u\ 

L 2 ~2 

is equal to A + if u < or — A_ if u > with \± as given in ( |5.12| ). 

From this, we have 1/2 Tr (T 2 + U 2 ) 1 ^ 2 = \u\L/A + T which, when substituting in equation (|1 . 7| ) 
(x2 in view of t), reproduces the ground state energy of the Example for Af > L if u > 
and Af < 2\A\ if u < (see equations ( p.20| ) and ( |5.23| )). Note that the equation ( |1.7| ) turns out to 




be equality in this example for < Af < 2L if u < and for \Af '/L — 1| < 1/2 if m > 0. 

Worth of mentioning here is the monotone behavior of E^> as a function of M. Monotonicity 
has been holden generically in all examples considered (satisfying bipartition or not) with the only 
exception being the Example pT3|.|I| for u > and < Af < L/2 (see equation ( p,19| )). This 



exception, as we believe, is related to the fact that the number of nuclei is related to the size of the 
sub-lattice: |^4| = iV n in this example. 

For u > there is an interval of values of Af, extending up to the half-filled band point 
Af = L, where E^ remains constant (see equations ( |3.8| ), ( |5.19| ) and ( p.28| )). We believe that 



this behavior is generically true and is related to the existence of first kind gap at Af = L for 
u > 0. As Af = N n + N e varies in this interval the minimum of E^ is attained with the nuclei 
number N n kept fixed. At Af = L all eigenvalues /i<0(/i<l/L for mean field hopping matrix) 
have been filled with electrons. As a consequence, if there is a second kind gap in the spectrum 
of H = T + u diag({cr x }), as usually do happen for typical hopping matrices T, the ground state 
energy density given by (|4.1| ) develops a kink at p = Af/L = 1/2 (see equation (|4.4j) ) giving rise a 
discontinuity in the chemical potential. 
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We finally discuss why there is no first kind gap if u < for the mean field models considered 
here. The key equation ( |3.11|) shows p,\ <•••<// N n < u < and \ij > if j > N n . So, for any 
number of particles Af = N n + N e , the minimum of energy is attained when all negative eigenstates 
are filled with electrons: N n = N e . The sum of these eigenvalues thus gives a smooth bulk (order L) 
contribution to Equations <KWi), (E23|), and ( ^35|) , show that e(p) = lim E {1 \Lp)/L 

is a smooth function of p, < p < 1 implying absence of first kind gap. However, if sub-dominant 
terms of are taking into account in equations ( |3.20| ) and ( |5.33| ), ^im [E^ l \L p) — Le(p)^j is 

continuous but non-different iable function of p at p = 1/2. This fact is the scar of a gap sweeped 
out by the long range hopping. Opposed to the result by Kennedy-Lieb, the mean field Falicov- 
Kimball model presents a phase transition at u — 0. Observe that there is no contradiction with 
item D of the introduction since the elements of the mean field hopping matrix are non-uniform 
in L. 



A Matrices of the Form al + (31 



Proof of Lemma \2.3\ . Clearly, the collection of matrices C form a vector space. This collection also 
forms a commutative algebra with the product of two matrices C = al + (31 and D = 7 J + Si 
given by 

CD = DC = cry / + (a6 + fry + f35R) 1 . 



From equation ( |A.1| ), D is the inverse of C, D = C -1 , if 

07 = 1 and a5 + fiy + [35R 

hold. Solving these two equations we get the coefficients of / and 1 of the item 2. 

To prove item 3 we notice that matrices of the form al + (31 are circulant matrices 



C = circfci, c 2 , 



CR 



Cl 


c 2 


c 3 . . 


■ Cr 


CR 


Cl 


c 2 . . 


■ Cr-i 


CR-1 


CR 


Ci . . 


■ Cr-2 


c 2 


c 3 


c 4 . . 


Ci 


= c 3 




• = Cr 


= (3. 



with 

ci = a + (3 and c 2 = c 3 = • • • = Cr — (3 . 
Circulant matrices can always be diagonalized by Fourier matrices F (see e.g. [PIP 

A = F- l CF 

where A = diag{Ai, A 2 , . . . , Xr} with 

Xj = Ci + c 2 uj r ~ 1 H h c R L) R R ~ 1){ ' J ~ l) . 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
(A.5) 

(A.6) 
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Substituting QA.4J) in ([0]) we get 

Ai = a + P R and A 2 

in view of the fact that 

j-l . , 2(i-l) 



A fi = a 



1 + c4 _i + w^"^ + • • • + wJT^ - ^ = (A.7) 
for all j = 2, ... , R. 

For completeness, let us verify that ( [A. 5 ) is true for matrices of the form C = ail + (31. Here 
we need only to Fourier transform the matrix 1 . Writing F as the matrix [1 ur^ • ■ ■ Ur^-i] with 
the vectors defined in Lemma 2~3 as columns and using ( A.7 ) give 

F^ 1 1 F = diag(.R, 0, ... 0) 

from which the eigenvalues of C can be read. 

Finally, item 4 follows from item 3 and this completes the proof of the proposition. □ 
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